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ABSTRACT 

We present 1420 MHz polarization images of a 2.5° x 2.5° region around 
the planetary nebula (PN) Sh 2-216. The images are taken from the Canadian 
Galactic Plane Survey (CGPS). An arc of low polarized intensity (size 0.2° x 0.7°) 
appears prominently in the north-east portion of the visible disk of Sh 2-216, co- 
incident with the optically identified interaction region between the PN and the 
interstellar medium (ISM). The arc contains structural variations down to the 
~1' resolution limit in both polarized intensity and polarization angle. Several 
polarization-angle "knots" appear along the arc. By comparison of the polariza- 
tion angles at the centers of the knots and the mean polarization angle outside 
Sh 2-216, we estimate the RM through the knots to be —43 ± 10 radm -2 . Using 
this estimate for the RM and an estimate of the electron density in the shell 
of Sh 2-216, we derive a line-of-sight magnetic field in the interaction region of 
5.0 ± 2.0 [iQ. We believe it more likely the observed magnetic field is interstel- 
lar than stellar, though we cannot completely dismiss the latter possibility. We 
interpret our observations via a simple model which describes the ISM magnetic 
field around Sh 2-216, and comment on the potential use of old PNe as probes of 
the magnetized ISM. 

Subject headings: planetary nebulae: individual (Sh 2-216) — ISM: structure - 
polarization — radio continuum: ISM 
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INTRODUCTION 



The diffuse Galactic synchrotron radiation provides a continuous background of radio 
emission which is intrinsically highly (up to ~70%) linearly polarized. This radiation is 
Faraday-rotated from the point of emission as it propagates through warm ionized gas in- 
terwoven with magnetic fields in the disk of the Galaxy; i.e., the angle 9 of the polarized 
component of the emission is rotated at wavelength A [m] by 



A9 = RMX 2 [rad], 



(1) 



where RM is the rotation measure [radm -2 ] and depends on the line-of-sight component 
of the magnetic field, B» \pG], the thermal electron density, n e [cm -3 ], and the path length, 
dl [pc], as 

RM = 0.81 J B\\n e dl [radm" 2 ]. (2) 
High-resolution radio polarization images at frequencies <3 GHz reveal the turbu l ent imprint 



of Fa r aday rotation on the diffuse polarized emission fe.g.. IWieringa et al 



199 



9: 



Gaensler et al. 



Haverkorn et al. 



2001 



2006 



b: 



Uvamker et a 



Schnitzeler et al 



1993 



Gray et al. 



2003; lHaverkorn. Katgert. fc de Bruyn 



2003a 



20071 ). The turbulent nature of the imprint is the 



product of the random component of the Galactic magnetic field and irregular electron- 
density distributions in the general interst ellar medium (ISM). Detailed studies and model- 



ing of the diffuse Galactic emission (e.g., ISpoelstralll984l ; lHaverkorn. Katgert. & de Bruvn 



2004a ) as well as statistical analyses of the RMs of polarized extragalactic sources (lHaverkorn et al 
2006al ) indicate that the scale size, or "cell" size, for variations in the magnetized ISM range 
from ~15 pc to 100 pc. Depth depolarization then results from the averaging of nonparallel 
polarization vectors from emission at different cells along the line-of-sight. Smaller-scale 
variations are also apparent in polarization images. Depolarization filaments, or "canals," 
with a width corresponding to one beam of th e observing instrument, indicate the pres- 



ence of very sharp gradients in RM (s ee, e.g., iGaensler et al.l l200ll ; lUyamker et al.l 12003 
Haverkorn. Katgert. &: de Bruynl l2004bl ) . If the gradient is so steep across the beam as to 
cause differential rotation of the polarization angle of ~90°, then complete depolarization 
occurs. Beam depolarization in images produced by aperture synthesis telescopes (typically 
with arcminute resolution) suggests a scale length for RM structures in the magnetized ISM 
of less than 1 pc. 

Interpreting the structu re seen in radio polariz ation images of the Galactic plane is 
largely left to modeling (e.g.. lHaverkorn et al.ll2004al ). as there is generally little correlation 
between the polarization structures seen in these images and the emission structures seen in 
total intensity images. Nevertheless, objects of known distance can be used to estimate the 
line-of-sight distribution of the magnetized ISM, revealing, for example, whether observed 



- 3- 



polarization structures are generated behind the known object, in the region between the 
object and the Sun, or perhaps i n the immediate vicinity of the object (see iGray et al. 
19991 ; lUyamker &: Landeckerl |2002| ) . Moreover, if small-scale depolarization structures can 



be isolated to a known object, then it may be possible to probe directly the properties of 
the magnetized ISM within the structures. 

H II regions and supernova remnants (SNRs) are the two most prevalent (discrete) con- 
stituents of the ISM as seen at radio wavelengths, and each class of object is detected in 
radi o polarization images up to a limiting distance determined by the "polarization horizon" 



( see lKothes fc Landeckerl 120041 ) . H II regions are detected by way of their depolarizing effects 
on background diffuse emission, while SNRs are simultaneously a source of polarized syn- 
chrotron emission and a Faraday "screen" which depolarizes background emission. However, 
neither class of object is a particularly good probe of the magnetized ISM. H II regions have 
very high electron densities and turbulent motions which produce tangled magnetic fields, 
a combination which leads to virtually complete beam depolarization across the region. For 
SNRs, complex models are needed to describe the physical parameters in the shock front at 
the interface between the rapidly expanding SNR and the ISM. 

Another class of object which may potentially serve as a better probe of the mag- 
netized ISM is planetary nebulae (PNe). Young PNe ar e relatively strong thermal radio 
emitters, resulting from high electron densities (see, e.g., iBains et al.ll2003l ). but they are 
also very compact (<Cl pc) and not yet interacting with the ISM. On the other hand, the 
shells of many old PNe are observed at optical wavelengths to interact with the ISM (see 
Tweedy fc Kwitterl Il996l ). Moreover, the ISM magnetic field appears to play a significant 
role in shaping the shells, as evidenced primarily by the visible "s t riping" or filamentar y 
structure of shell gases (see lTweedy. Martos. fc Noriega-Crespolll995l ; ISoker fc Zuckerlll997l ). 
Theoretical treatments show that interactions between PNe and the ISM are an important 
consideration in the evolution of PN systems moving at even modest speeds (>5 kms~* ) 
with respect to the ISM (see, e.g.. lSoker &: Dganilll997l ; IWareing. Zijlstra. fc O'Brien! 120071 ). 
If the conditions are right in the interaction region between the PN and the ISM, we may 
expect to see the Faraday signature of old PNe in radio polarization images. Such a signa- 
ture has been id entified fo r the n earby PN Sharpless 2-216 (Sh 2-216) and was first described 
by one of us in lUyamkerl (120041 ) . In this paper, we describe in detail the Faraday-rotation 
structure in the shell of Sh 2-216. 

We present radio polarization images at 1420 MHz of a 2.5° x 2.5° region of the Galactic 
plane around the position of Sh 2-216. The images are taken from the Canadian Galactic 
Plane Survey (CGPS). In § [2j we summarize the pertinent properties of Sh 2-216. In § [3J we 
describe briefly the preparation of the images. In § [U we describe the structures observed on 
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the visible disk of Sh 2-216 in both polarized intensity and polarization angle, and estimate 
the RM through the shell of Sh 2-216. In § El we derive the magnetic field in the shell of 
Sh 2-216, and interpret the structures and RMs in the context of an interaction between the 
PN and the ISM. We also discuss the possibility that the observed structures are produced 
by the stellar field of the host white dwarf or its progenitor. Finally, in § |6l we summarize 
our conclusions. 



2. THE PLANETARY NEBULA SH 2-216 



Sh 2-216 is the closest known PN. At a distance of 129 pc (Harris et all 120071 ). its 1.7° 
angular diameter translates to a physical diameter of 3.8 pc, making it also one of the 
largest and oldest PNe. The most conspicuous feature of Sh 2-216 at optical wavelengths is 
its bright eastern rim, denoting an interaction between the expanding and moving PN and 
the ISM. The location of the interaction region^ appears to be consistent with the observed 
displacement of the host white dwarf from the center of the PN; i.e., the enhanced emission in 
this region is a consequence of the additive velocity of the nebular expansion in all directions 
and the underlying eastward moti on of the PN sy stem relative to the ISM. The expansion 
velocity is very low (< 4 kms -1 ; iReynoldd 119851 ). indicating that the ISM pressure (with 
dynamic, magnetic and cosmic-ray components) is nearly equal to the PN ram pressure. The 
velocity (o n the plane of the sk y) of the system relative to the ISM is estimated also to be 
~4 kms -1 (ITweedy et al.lll995l ). 



The thin filamentary structures observed in Ha in the interaction region, together with 
the more subtle, and wider, filamentary structure observed in N II across the face of the 
PN, qualitatively su ggest that the ISM magnetic field is shaping the morphology of Sh 2-216 
( ITweedy et al.lll995l ). Using estimates for the elec tron density w i thin t he PN (n e ~ 5 cm -3 ) 
and the mean ISM magnetic field (B ~ 5 /xG), ITweedy et al.l ( 119951 ) show that the ISM 
magnetic pressure is about twice that of the dynamic pressure, and thus likely a dominant 
factor in the shaping. The strength of the PN magnetic field is expected to be negligible at 
large rad ii (r ~ 1 pc) from the host white dwarf, assum ing the field decreases as B ~ r~ 2 



v see, e.g. 



within filaments due to compression (see 



Soker 



Vlemmings. Diamond, fc van Langeveldeil2002h. but may be amplified significantly 



2002 



Huggins fc Manley 



20051 ). We discuss 

the strength and orientation of the magnetic field in the outermost regions of Sh 2-216 in §0 



1 We refer to the bright eastern rim as "the interaction region" throughout the paper, though other parts 
of the shell of Sh 2-216 may also be interacting with the ISM. 
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OBSERVATIONS AND IMAGE PREPARATION 



The radio polarization d ata presented in th is paper were obtained at 1420 MHz (A = 
21cm) as part of the CGPS (Taylor et al.ll2003l ) using the synthesis telescope (ST) at the 
Dominion Radio Astro physical Observatory (DRAO). The ST is described in detail by 
Landecker et al.l (120001 ) . Images are produced in each CGPS field for the two orthogonal 
linear polarization states, Stokes-Q (Q) and Stokes-U (U), as well as Stokes-I (total inten- 
sity), from data in each of four 7.5 MHz continuum bands centered on 1406.65, 1414.15, 
1426.65 and 1434.15 MHz, respectively. Images in Stokes- V, nominally representing cir- 
cular polarization, are presently dominated by instrumental errors, and are of significantly 
lower value. The DRAO ST is sensitive at 1420 MHz to emission from structures with an- 
gular sizes of ~1° (corresponding to the shortest, 12.9 m, baseline of the ST) down to the 
resolution limit of ~1' (corresponding to the longest, 61 7.1 m, baseline). I n total intensity, 
data from the Effelsberg 21-cm Radio Continuum Survey (IReich et al.ll 199 71 ) are added to the 



band -averaged ST data to provide information on the largest spatial scales (see ITaylor et al. 



20031 ). In Q and U, data from t wo single-antenna surv eys of the northern sky at ~1.4 GHz, 
namely the DRAO-26m survey (jWolleben et al.l 120061 ) a nd the Effelsberg Med ium Latitude 
Survey, are added to the band-averaged ST data (see lLandecker et al.l 120081 ). All CGPS 
images presented in this paper were produced using band-averaged data. 



CGPS data calibration and processing procedures are described in detail in ITaylor et al. 



(120031 ). Here we summarize for the reader the general practice, emphasizing procedures re- 
lated specifically to polarization. The complex antenna gains for the pointing centers in 
each CGPS field are calibrated by observing a compact calibrator source, either 3C 147 or 
3C 295, at the start and end of each observing session. The polarization angle is calibrated 
using the polarization calibrator source 3C 286. Amplitude and phase variations encountered 
during individual observing sessions (on time scales down to ~2 hr) are determined during 
the processing of the total intensity data, and are applied also to the Q and U data. Addi- 
tional processing, needed to remove the effects of strong sources both inside and outside the 
primary beam of t he ST antenn as, is accomplished using routines developed especially for 
the DRAO ST (see IWillid 119991 ). The instrumental polarization, which is corrected on-axis 



in the sequence above, varies across the primary beam of the ST antennas due to cross- 
polarizatio n of the (nomin ally orthogonal) receiver feeds and the effects of the feed support 



struts (see iNg et al.l 120051 ) . The result is "leakage" of unpolarized radiation, seen in total 
intensity, into Q and U. We have employed two different methods at DRAO to correct for 



2 Note that the frequency corresponding to the midpoint of the four continuum bands is 1420.4 MHz, the 
neutral hydrogen spin-flip frequency. The 5.0 MHz band about this frequency is allocated to the 256-channel 
spectrometer (see lTavlor et al.ll2003l ) 
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the wide-field instrumental polarization. Both methods were used to calibrate the various 
CGPS fields appearing to some degree in the 2.5° x 2.5° region presented in this paper. In 
the first method, we derived the "average" leakage pattern in Q and U across the primary 
beam of the ST antennas, and subtracted from each of the Q and U imag es the "leakage 
image" for total intensity into Q and U, respectively (see lTaylor et al.ll2003l ). In the second 
(and newer) method, we derived the leakage patterns for each of the ST antennas separately, 
and subtracted the complex leak age pattern creat ed by each pair of antennas directly from 
the Q and U visibility data (see iReid et al.ll2008l ). The residual instrumental polarization 
error after the on-axis and wide-field calibration is similar for each method in each processed 
field, increasing from ~0.3% root-mean-square (rms) at the field pointing center to ~1% 
at the field edge (p = 75'). The rms error is reduced further in the mosaicing process. 
The newer wide-field correction significantly reduces artifacts in Q and U around bright 
total-intensity sources (i.e., with flux densities >100 mJy). No artifacts are seen above the 
estimated ~0.34 mJybeauT 1 (~0.086 K) noise level in the 2.5° x 2.5° region of the CGPS 
presented in this paper. 



4. RADIO POLARIZATION IMAGES OF PN Sh 2-216 

In Figured] we show images of the 2.5° x 2.5° region around PN Sh 2-216 in both optical 
intensity at R-band (A = 6570 nm) and total radio intensity at 1420 MHz (A = 21cm). The 
optical image is taken from the Digitized Sky Survey (DSS) and the radio image from the 
CGPS. The images are presented in Galactic coordinates and centered on the position of the 



host white dwarf LS V46°21 (I = 158.49°, b = +0.47°; see lKerber et al.ll2003f ). Note that the 
center of t he visible disk of S h 2-216 is offset «24' to the Galactic south-west of the white 
dwarf (see iTweedv et aDl995h . For the optical image, we adjusted the range of intensities to 



highlight extended emission. For the total radio intensity image, we removed point sources 
leaving only extended emission. There is a clear enhancement in both images across much 
of the face of Sh 2-216, relative to the surroundings, but the enhancement is most intense 
along the (Galactic) north-eastern rim; i.e., in the interaction region between Sh 2-216 and 
the ISM. Treating the enhancement in the radio image as thermal emission from the shell of 
Sh 2-216, we can estimate the thermal electron density in the interaction region of the PN 
(see § ED- 

In Figure [2] we show the polarized intensity (P = ^/Q 2 + U 2 — (1.2cr) 2 , where the last 
term gives explicitly the noise bias correction) and polarization angle (9p = h arctan U/Q) 
images at 1420 MHz for the 2.5° x 2.5° region in the CGPS around PN Sh 2-216. The 
polarization images contain several interesting features on a variety of angular scales. The 
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most notable feature is a low-polarized- intensity arc ~0.15° wide and ~0.7° in length, co- 
inciding with the north-east portion of the visible disk of Sh 2-216. The reduced intensity 
and distinct shape of the arc indicate that its appearance is due to the effects of Faraday 
rotation: specifically (1) localized beam depolarization within the arc of the background 
diffuse synchrotron emission, and/or (2) cancellation of the background emission, Faraday- 
rotated within the arc, by foreground synchrotron emission. For background/foreground 
cancellation to play a significant role, the polarized foreground emission must be a reason- 
able percentage of the total polarized emission in the direction of Sh 2-216. Galactic models 
for synchrotron emission predict for the 129 pc foreground t oward Sh 2-216 only ~0.06 K at 



1420 MHz (e.g.. iBeuermann. Kanbach. fc Berkhuijsenlll985l ). Even if this emission is highly 



(i.e., ~70%) polarized, we expect just ~0.04 K of polarized emission in the foreground; i.e., 
<10% of the 0.47 ± 0.06 K total polarized emission seen outside the visible disk of Sh 2-216 
(and at b > +0.5°). We conclude, therefore, that beam depolarization plays a larger role in 
reducing polarized emission over the arc than background/foreground cancellation. Struc- 
tural variations observed within the arc in both polarized intensity and polarization angle on 
angular scales down to the resolution limit (~1') further suggest that beam depolarization 
is responsible for the appearance of this feature. Since the length and location of the arc are 
very similar to the optically bright rim denoting the PN-ISM interaction region, it would 
seem there is a physical connection between the conditions and processes in this region which 
give rise to enhanced optical emission and those which lead to sharp gradients in RM. We 
discuss the RM structure within the arc in § 14.11 

Aside from the prominent north-east arc, do we see other signatures of Sh 2-216 in the 
polarization images? The circle representing the visible disk of Sh 2-216 in Figure [2] draws 
our attention to two suggestive details in the northern half (b > +0.5°) of the images: (1) 
the appearance of a second low-polarized-intensity "arc" ~0.2° wide and ~0.4° in length, 
located in the north-west portion of the visible disk of Sh 2-216; and (2) the increased range 
of polarization angles on the visible disk of Sh 2-216 compared to the surroundings. The 
small-scale structural variations in polarization angle within the north-west arc indicate that 
beam depolarization is responsible to at least a moderate degree for the reduced emission 
in this feature. If this second arc is indeed associated with Sh 2-216, then the conditions for 
sharp RM gradients in the shell of the PN may not be confined to the optically-identified 
interaction region. Moreover, a comparison of the range of polarization angles seen inside 
the visible disk of Sh 2-216 (-82° to +27°, rms » 17°) with those seen outside (-17° to 
+32°, rms ~ 6°) suggests that the conditions for moderate RMs are present throughout the 
shell of Sh 2-216. We present a simple model for the observed polarization structures on the 
visible disk of Sh 2-216 in § 



In contrast to the smaller-scale structures seen on the visible disk of Sh 2-216 in the 
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northern half of the images (b > +0.5°), the southern half of the images (b < +0.5°) is 
dominated by "bands" of relatively low polarized intensity, 0.1°-0.3° wide, which stretch 
approximately east-west across the region. The bands have no counterpart in total intensity 
The boundary between north and south is clearly marked in the polarization angle image by 
a jagged line over which the angle changes very rapidly. On close inspection of the polarized 
intensity image, this line corresponds to a narrow (~1') channel within the northernmost 
band of virtually zero polarized emission. Changes in the polarization angle across "cells" 
< 3' in size are seen, to differing degrees, throughout the bands. The bands appear to be 
part of a large-scale complex which depolarizes the background diffuse emission, most likely 
before it reaches the position of Sh 2-216. In the less likely scenario that the complex sits 
between the Sun and Sh 2-216, any polarization signature imprinted on the background by 
Sh 2-216 is lost. In either case, the positioning of the bands on the sky south of the north-east 
arc associated with Sh 2-216, and other apparent features on the northernmost portion of 
the visible disk of the PN, would seem to be fortuitous. 



4.1. Rotation-Measure Structure in the North-East Arc 



The polarization angle of the relatively bright emission outside the visible disk of 
Sh 2-216 (and at b > +0.5°) has a mean value of +7° and rms variations of only 6°. Along 
the prominent north-east arc, the polarization angle is observed to change rapidly across the 
perimeters of roughly elliptical "knots." The angle inside the perimeters changes more slowly 
and, indeed, plateaus at the centers of the knots. In Figure [3] we show a small 0.4° x 0.4° 
region in polarization angle around the north-east arc and identify eight discrete knots. We 
define as the center of each knot the position of the pixel showing the maximum clockwise 
(see below) deviation from the background (+7° ± 6°) value. In Tabled] we give the mean 
value of the polarization angle in each knot. The mean was estimated over an area corre- 
sponding to the area of the resolving beam (10 pixels, see Figure [3]), excluding, in the cases 
of knots 2 and 4, pixels which differed from the 10-pixel mean by more than 2a. Table [T] 
shows that the polarization angle of the emission emerging from the knots is rotated signif- 
icantly with respect to the background emission. The weighted mean polarization angle at 
the centers of the knots is +78° ± 22°. If we assume that emission with polarization angle 
+7° ± 6° is incident on the far side of each knot, and for the moment ignore foreground 
emission, then the incident emission is Faraday rotated in the knots by A6 = —109° ± 23°. 
We infer negative, i.e., clockwise, rotation by tracing polarization angles from the outside 
edge of the arc to the center of any knot. The trace shows that the polarization angle (first) 
decreases through negatives values. At five of the eight knot perimeters, the polarization 
angle jumps from —90° to +90°, and then continues to decrease to its center value. Since 
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foreground emission probably cannot be ignored at the ~10% level, we must estimate the 
maximum deviation expected in the observed polarization angle if, by chance, the foreground 
emission is rotated 45° relative to the emission emerging from the knots. (Note that fore- 
ground emission rotated 90° relative to the background leads to a maximum reduction in 
polarized intensity, but no net rotation in polarization angle.) Assuming complete beam de- 
polarization at the perimeters of at least some of the knots, in particular knots 1 and 4, we 
estimate the polarized foreground emission to be 0.046 ±0.012 K, consistent with the values 
predicted by Galactic synchrotron models. Using 0.177 ± 0.044 K for the mean observed 
(i.e., emerging plus foreground) polarized emission at the centers of the knots (see Tabled]), 
we estimate a maximum deviation of 7°. Adding this in quadrature to the 23° statistical 
uncertainty gives a standard error for the measured rotation through the knots of 24°. For 
a center wavelength of 21.12 cm (see § EJ), a rotation of —109° ± 24° gives (via Equation 1) 
RM = -43 ± 10 radm- 2 . 

The emission emerging from the centers of the knots, where the polarization angles 
plateau, is likely higher than the 0.177 ± 0.044 K value given above, but polarization-angle 
variations over the resolving beam lead to reduced polarized intensity even inside the knot 
perimeters. While these variations are reflected in the range of angles observed in each knot 
(see Table [1]), we nevertheless believe our estimated mean RM reflects a real systematic 
rotation of the background emission as it passes through the knots. The consistent clockwise 
rotation of the polarization angle observed in moving from outside the north-east arc toward 
the center of any knot strengthens this assertion. 

We tried to estimate the RM through the knots using the four-band data from the ST 
but failed. Given the typical uncertainty in polarization angle in the band-averaged image, 
and noting that RM « —43 radm~ 2 gives a difference in rotation angle of only ~5° over 
27.5 MHz (AA = 0.41 cm), the failure is not surprising. 



We can derive the line-of-sight component of the magnetic field through the knots in 
the north-east arc of Sh 2-216 using —43 ± 10 radm" 2 as an estimate of the RM in the 
knots, and using estimates of the thermal electron density in and path length through the 
interaction region (see Equation 2). The electron density in Sh 2-216 can be calculated 
from emission measure (EM = J n e 2 dl) determinations, made independently at optical 
and radio wavele ngths. Ba s ed on their measured Ha intensity and gas temperature (T e = 



Sh 2-216 of EM « 42 cm 6 pc. Using a brightness temperature of Tj, = 0.11 ± 0.02 K for the 



5. 



DISCUSSION 




value for the emission measure over 
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thermal radio emission in the interaction region of Sh 2-216 (obtained via comparison of on- 
source and off-source temperatures in Figure and the same gas temperature, we estimate 
a value for the interaction region of EM = 69±13 cm -6 pc. If we assume for the moment that 
electrons are uniformly distributed over the approximately spherical volume of Sh 2-216, and 
use an average path length through the sphere of Al = f i?pN ~ 2.5 pc, then the optically- 
determined EM gives a mean electron density over Sh 2-216 of n„ ~ 4.1 cm -3 . In some 
contrast, the radio-intensity-determined EM gives, for a path lengths through the interaction 
region of Al = 1.1 ± 0.3 pc, n e = 7.9 ± 1.3 cm -3 . The factor ~2 increase in the electron 
density in the interaction region compared to the mean value over the entire PN is reasonable, 
since material is stacking up in the interaction region (see iTweedy et al.lll995l ). However, 
n e = 7.9 ± 1.3 cm" 3 still represents a mean over the interaction region. The Ha images of 



Tweedy et al.l (119951 ) show small-scale filamentary structures in the interaction region with 



localized factor 1.5-2 enhancements in EM relative to the mean. The filaments correspond 
approximately in both location and size to the polarization-angle knots. Since we cannot 
confirm the physical association between the filaments and the knots, we conservatively 
assume an E'M-enhancement of 1.5 ± 0.5, and estimate the electron density in the knots 
to be n e = 9.7 ± 2.3 cm" 3 . Using RM = -43 ± 10 radm" 2 , n e = 9.7 ± 2.3 cm" 3 and 
Al = 1.1 =h0.3 pc, we derive a line-of-sight magnetic field through the knots in the interaction 
region of B\\ = 5.0 ± 2.0 /xG. Since the RM is negative, this field is directed into the plane 
of the sky. 



5.1. An ISM Origin for the Magnetic Field in the Shell of Sh 2-216 



Is a ~5 /iG line-of-sight magnetic field reasonable for the ISM around Sh 2-216? Since 
there is no direct measurement of the ISM magnetic field around Sh 2-216, we estimate the 
local field from what is known generally about the Galactic magnetic field. T he Galactic 
magn etic field is concentrated in the disk and has two components (see, e.g., iBeck et al. 
19961 1: a large-scale or regular component (B reg ), which follows the spiral arms, and a 



small-scale or random component (B ran ). B reg in the local spiral arm is found, using po- 
larized radio sources and the polarization of starlight, t o be directed toward I ~ 85° (e.g., 
Rand fc Lyndll994l ; iHeilea Il996bl : iBrown fc Taylor! 120011 ); i.e., clockwise as viewed from the 
Galactic north pole. Fluctuation cell sizes fo r B ran are estimated to be 50-100 pc (see 
Rand fc Kulkarnilll989l ; lOhno fc Shibatal 119931 ). The ratio of the strengths of the random 
and regular components of the Galactic field, B ran /B reg , can be obtained directly from 



3 The path length Al = 1.1 ± 0.3 pc corresponds to the mean of the line-of-sight chord lengths through a 
1.9-pc radius sphere at the positions of the eight knots. 
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starlight po larization data and synchrotron pola rization data using t he model presented in 
Burn! (119661 ). The starlight polarization data of iFosalba et al.l (120021 ) give for a large sam- 
pl e of stars covering a ll Gala ctic longitudes B ran / B reg 1.3. Using stars from the sample 
of iMathewson fc Fordl (119701) in the range 120° < / < 180°, and a modified version of the 
Burn model, Heilesl ( 1996a ) finds B ran /B reg 1.5. For synchrotron emission just north 
of the visible disk of Sh 2-216 (7 = 158.5°), we measure a fractiona l linear polarization of 
p = 0.27 ± 0.03, close to the maximum value found by ISpoelstral (119841 ) for the diffuse 
emission in the Galactic plane. With a spectral index a = —0.44 ± 0.04 (S oc v a ) between 
408 MHz and 1420 MHz for the synchrotron emission in the CGPS region around Sh 2-216, 
we get for the intrinsic value of the fractional linear polarization p max = 0.68 ± 0.01 (see 



Ginzburg &: Svrovatskiilll965l ). and thus obtain B ran /B reg = 1.51 ±0.13, consistent with the 



Heilesl (Il996al ) starlight estimate. Using a value B tot = 4.2 jiG fo r the average azimuthal 
field strength (B^ ot = B 2 reg + B^ an ) in the local arm (IHeilesI Il996al ) and B ran /B reg = 1.51, 
we estimate B reg fh 2.3 fiG and B ran fh 3.5 fiG. The maximum magnetic field at any point 
in the local arm is then achieved if, by chance alignment, the random field lies parallel to 



the regular field; i.e., B n 



Breg ~t~ B T 



5.8 /xG. At the longitude of Sh 2-216, both 



the average field (B tot = 4.2 /iG) and maximum possible field (B max = 5.8 /zG) lie largely 
in the plane of the sky, and run from Galactic east to west. The maximum field along the 



line-of-sight, where B\\ reg 
plane of the sky. 



0.6 /zG, is B\\ 



B\ 



rey 



+ B r 



4.1 fiG, directed into the 



In light of this brief overview, we conclude that an intrinsic ~5 fiG line-of-sight magnetic 
field in the ISM at the position of Sh 2-216 is unlikely. Nevertheless, our observations can be 
used to comment further on the structure of a proposed ISM field in the interaction region 
as well as other locations in the shell of Sh 2-216. Our estimate of the line-of-sight magnetic 
field in the interaction region is based on the maximum RM as seen through knots in our 
polarization angle image. The RMs outside the knots are apparently much lower. The 
sharp RM gradients over the knot perimeters must be the result of either a rapid change 
in the electron density or the line-of-sight magnetic field, or both. As we previously noted, 
the polarization-angle knots appear to be associated with narrow Ha filaments observed in 
the interaction region. The sharp edges of the filaments, which denote a rapid change in 
EM (and thus electron density), naturally explain RM gradients across knot perimeters. 
The magnetic field need not change significantly across the interaction region. Realistically, 
however, the magnetic field is probably affected by turbulence in the hot gas (see § 15.1.11) . 

If we look west of the north-east arc, toward the center of the visible disk of Sh 2-216, 
we continue to see polarization angles significantly different from the +7° ± 6° observed 
outside the PN (see Figure©). Though we don't see prominent structures in this "interior" 
region, we do see some localized polarization-angle structures. These structures are roughly 
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coincident with low-level enhancements in Ha and N II (see iTweedy et al.l Il995l ) and total 
radio intensity (see Figure [1^). Localized electron-density enhancements may therefore be 
responsible for both the knots in the interaction region and the more extended structures seen 
across the western portion of the face of Sh 2-216. Indeed, these two apparently different 
structures may arise from similar underlying structures, seen e dge-on in the case of the 



knots, and face-on in the case of the extended structures (see, e.g.. lTweedy et al.lll995l ). The 
increased path length through the shell in the interaction region would explain, at least in 
part, why the EMs and RMs in the filaments/knots are larger than those in the extended 
structures across the face. A decrease in the line-of-sight component of the magnetic field, 
moving west from the north-east edge of Sh 2-216 toward the center of the visible disk, could 
also account for some of the difference (see § 15.1.11) . 

A second low-polarized-intensity arc appears at the north-west edge of the visible disk of 
Sh 2-216 (see Figure^). Small-scale variations in polarization angle within this arc indicate, 
as in the north-east arc, the presence of sharp RM gradients. However, unlike the north-east 
arc, there are no plateaus (i.e., multi-pixel regions of roughly constant polarization angle) over 
which we can confidently estimate some deviation from the outside +7° ±6°. Consequently, 
we have no means of estimating the magnitude of the RM through this arc. Nevertheless, 
there is some indication of the sign of the RM. Moving south from the edge of the north-west 
arc, the polarization angle (on average) increases, implying positive RMs. To substantiate 
th is finding, we broke the no rth-west arc into three north-south slices, and used the approach 



of IWolleben &: Reich! (120041 ) to estimate RM together with three other parameters (degree 
of depolarization, foreground polarized intensity and foreground polarization angle). We 
found positive RMs for each slice, even when we varied the other parameters away from 
their "best-fit" values. The positive RMs indicate that the magnetic field in the north-west 
arc is directed out of the plane of the sky. If the north-west arc is associated with Sh 2-216, 
and the ISM field is responsible for the observed RMs in both the north-east and north-west 
arcs, then the intrinsic field must be deflected around the PN. 



5.1.1. A simple model for the ISM magnetic field around Sh 2-216 

For the ISM magnetic field to simultaneously account for the negative RMs observed in 
the north-east arc and the positive RMs observed in the north-west arc, the intrinsic field 
must bend significantly around the shell of Sh 2-216 such that it has a ~5 /iG line-of-sight 
component into the sky on the east edge of the PN and a non-zero line-of-sight component 
out of the sky on the west edge. This is exactly what we might expect in the following 
scenario (see Figure H]): The intrinsic ISM field around Sh 2-216 is described by the 4.2 /iG 
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azimuthal component of the Galactic magnetic field (see lHeileall996a[ ). which, at / = 158.5°, 
runs in the local arm from Gala ctic east to west and intersects the plane of the sky at 
15.5° ± 4° (IBrown &: Taylor! 1200 ll ). The intrinsic field is compressed and deflected by the 
expanding and movin g PN, since it can diffuse only slowly into the partially ionized shell (see 
Soker fc Dganilll997l ). The three-dimensional motion of Sh 2-21 6 is fully characterized by the 



(Galactic) north-we st-directed motion of the host white dwarf (ICudworth fc Reynoldslll985 



Tweedy et al.lll995l ) and a line-of-sight motion into the plane of the sky (see below). The 
motion in the plane of the sky is of less importance for our observations than the motion 
along the line-of-sight, though the full three-dimensional picture is important for interpreting 
the alignmen t of the filaments in t he interaction region and wider structures across the face of 
Sh 2-216 (see lTweedy et al.lll995l ). At the far side of Sh 2-216, the line-of-sight motion drags 



the intrinsic field away from the observer. The result is a deflected field around Sh 2-216 
which has on the east edge of the PN a line-of-sight component directed into the sky and 
on the west edge a line-of-sight component directed out of the sky. In the center portion 
of the PN, the field lies largely in the plane of the sky. The line-of-sight component of the 
field on the east edge is slightly larger in strength than the intrinsic field itself, while that 
on the west side is lower. We estimate the l ine-of-sight motion of Sh 2-216 by consulting 
the Wisconsin Ha Mapper (WHAM) survey JHaffner et alll2003h . The WHAM data show 
that the Ha emission from Sh 2-216 peaks at velocity +5 ± 1 kms -1 relative to the local 
standard of rest. (Note that the sign for velocity is opposite that of RM; i.e., a positive 
velocity signifies motion into the plane of sky while a positive RM signifies a magnetic field 
out of the plane of the sky). The emission surrounding Sh 2-216 peaks at velocities near 
zero, indicating that the +5 ± 1 kms -1 is indeed relative to the surrounding ISM. 

The magnetic field described by our model represents only the "smooth" component of 
the deflected ISM field. The magnetic field in the shell of the PN will also have a turbulent 
component due to motions in the hot gas. The turbulent component contributes in part 
to the beam depolarization we observe in the interaction region. The smooth component is 
responsible for the systematic RM observed through the knots. 



5.1.2. Are old PNe good potential probes of the magnetized ISM? 

We have asked in this section whether or not the magnetic field derived from the observed 
RMs is reasonable for the ISM around Sh 2-216. If we were, instead, to concede that the 
derived field is native to the ISM, then we could ask the question: Are old PNe, such 
as Sh 2-216, good probes of the intrinsic ISM field? The simple qualitative model presented 
above for Sh 2-216 suggests that we can learn something about the intrinsic field surrounding 
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this PN. With a more comprehensive (three-dimensional) model of the PN-ISM interaction, 
it may be possible to work out quite accurately the strength and orientation of the intrinsic 
field . Since a large percentag e of old PNe show PN-ISM interactions similar to Sh 2-216 



see 



Tweedy fc K witter! Il996l ). it is perhaps reasonable to assume that the conditions for 
detectable Faraday rotation are present in the shells of many old PNe in the nearby Galaxy. 
Given a good model of the interaction in each case, it should then be possible to determine 
the intrinsic field at many locations. 

There are two significant drawbacks to consider before we declare old PNe good potential 
probes of the magnetized ISM: (1) It may not be possible in many cases to construct a 
good model of the PN-ISM interaction, due either to large uncertainties in the physical 
parameters (e.g., ISM and PN particle densities, space velocity of the PN) or to the overall 
complexity of the interaction. The qualitative model presented above for Sh 2-216 does not 
comment on either the degree to which the intrinsic field around the PN is compressed or 
the maximum angle with which the intrinsic field is deflected. Both of these quantities are 
necessary in order to more accurately determine the intrinsic field from the line-of-sight 
field. Three-dimensional magneto-hydrodynamic (MHD) simulations can perhaps be used 
to demonstrate the sensitivity of the intrinsic-field determination to various parameters. (2) 
The Faraday signature of the PN is at the mercy of fluctuations in the warm ionized ISM 
as well as turbulent structures (e.g., H II regions) which may lie along the line-of-sight. 
The dark bands that run through the approximate midpoint of Sh 2-216 (almost) completely 
depolarize the diffuse background emission. If Sh 2-216 were located ~0.5° south of its actual 
position, then the distinct signature of the north-east arc would be destroyed. Given the 
prevalence of turbulence in the ISM, this point is a significant concern. In fact, of the six PNe 
in the CGPS region known to interact with the ISM, only two, including Sh 2-216, are seen 
in polarization. (The other, namely DeHt5, is the subject of a subsequent paper.) Targeted 
polarimetric observations of old PNe at sub-arcminute resolution and at multiple frequencies 
in the range 1-3 GHz are necessary to better establish the potential of these objects as good 
probes of the magnetized ISM. 



5.2. A Stellar Origin for the Magnetic Field in the Shell of Sh 2-216 



Can a ~5 fiG magnetic field in the interaction region be attributed to the host white 
dwarf near the center of Sh 2-216? The magnetic fields of white dwarfs have been measured 
via spectropolarimetric observations of optical absorption lines, but only recently have the ob- 
servations had the sensitivity to detect kilog auss fields. The studies thus f ar have focused on 
either fully evolved (compact) white dwarfs (lAznar Cuadrado et al.ll2004j ) or central stars of 
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relati vely young PNe which are still transiti oning; to white dwarfs (Jordan. Werner, fc O'Took 
2005 ). In the case of evolved white dwarfs, Aznar Cuadrado et ali ( 2004 ) found for a sample 
of 12 stars only three which had d etectable magnetic fi elds in the range 2-4 kG, a detection 
rate of 25%. On the other hand, iJordan et ajj (120051 ) found for each of a selection of four 
transition stars magnetic fields of 1-3 kG, a detection rate of 100%. Though the number 
statistics for both cases are relatively poor, this pair of results suggests that the magnetic 
fields of transition stars are present not in their degenerate cores but rather their extended 
envelopes, since magnetic flux is apparently lost dur ing white dwarf evolution (i.e., during 
collapse from stellar radii in the IJordan et al.l 120051 sample of 0.14-0.3 Rq to white dwarf 
radii of ~0.012i? Q ). The magnetic fields of the central stars of old PNe have not been mea- 
sured. Given their intermediate radius, e.g., 0.05 R & fo r the "central" star in Sh 2-216 (based 
on the luminosity and effective temperature given in iRauch et al.l 120071 ) . measurements of 
the magnetic fields of the central stars of old PNe could lead to an improved understand- 
ing of field evolution in white dwarfs. We point out for completeness that a small fraction 
(^1 0%) of white dwarf s are observed to have magnetic fields at the 1 MG level or higher 



sec 



Liebert et al.l 120031 ) . but these stars tend to have masses (~O.9M ) much higher than 
typical white dwarf masses (0.48-0 . 65 MJ, and may come fr om magnetized progenitors such 
as peculiar (Ap) stars (see lLiebertlll988l : iLiebert et al.l 120031 ) . 



For present purposes, we assume the magnetic field in the envelope around the contract- 
ing "central" star in Sh 2-216 to be accurately represented by the magnetic field {B^ z — 
1.8 kG ) measured at the radii (r avg = 0.21 R Q ) of the central stars in the IJordan et al. 



(120051 ) sample. With some knowledge of the large-scale magnetic field geometry, we can 
then estimate the field at large radii; namely, in the shell of Sh 2-216. Unfortunately, at 
this time, neither observations nor theory form a complete picture of the magnetic fields 
in PNe. Magnetic field measurements of maser spots in precur sor (AGB) circumstella r 
envelopes suggest a radial dependence of the field (B ~ r~ 2 ; see IVlemmings et al.ll2002l ). 



~ r 3 ; see 



Vlemmines. van Laneevelde. & Diamond 


2005 V 


structures observed by 


Hue:e:ins & Manlev 


(2005) 


Vlemmines, Diamond. & Imai 


( 


2006) of the mag 



In contrast, the geometry of filamentary 



from AGB star W43A, sugge st the dominance of toro idal fields (B ~ r" 1 ), consistent with 
the theoretical framework of I Chevalier fc Luol (I1994J ). If either radial or poloidal geome- 
tries hold for Sh 2-216, then the magnetic field in the interaction region (»1.0 pc from the 
host white dwarf) would fall well below our i?M-estimated ~5 /iG line-of-sight field. On the 
other hand, if a toroidal field holds, then the magnetic field in the interaction region could be 
~8 /iG. Given the spherical symmetry of the shell of Sh 2-216, a large-scale toroidal magnetic 
field for this PN, invoked generally to explain non-spherical (e.g., bipolar, elliptical) symme- 
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tries in young PNe, is unlikely. However, localized enhancem ents of the i nternal magnetic 
field, due to compression in dense knots or filaments (see, e.g., ISokerl 120021 ; ISoker fc Kastner 
20031 ). are possible. Thus we cannot completely dismiss the possibility of a ~5 fiG internal 



field in the shell of Sh 2-216. Detailed MHD simulations need to be done in order to better 
understand the magnetic field geometry in PNe. 



6. CONCLUSIONS 

Here we give a summary of our results and conclusions: 

1. We presented 1420 MHz polarization images for the 2.5° x 2.5° region in the CGPS 
around the PN Sh 2-216. 

2. A low-polarized-intensity arc, 0.2° x 0.7° in size, appears in the north-east portion 
of the visible disk of Sh 2-216. The arc is coincident with the optically-identified interaction 
region between the PN and the ISM. 

3. A second low-polarized-intensity arc appears in the north-west portion of the visible 
disk of Sh 2-216. 

4. The north-east arc contains structural variations down to the resolution limit in 
both polarized intensity and polarization angle. Several polarization-angle "knots" appear 
along the arc. 

5. Via comparison of the polarization angles at the centers of the knots in the north-east 
arc and the mean polarization angle outside Sh 2-216 (and above b ~ +0.5°), we estimated 
the RM through the knots to be —43 ± 10 radm -2 . 

5. Using this estimate for the RM and an estimate of the electron density in the shell of 
Sh 2-216, we derived a line-of-sight magnetic field in the interaction region of 5.0 ± 2.0 /zG. 

6. We believe it more likely the derived magnetic field is interstellar than stellar, though 
we cannot completely dismiss the latter possibility. We interpret our observations via a simple 
model which qualitatively describes the ISM magnetic field around Sh 2-216. 

7. It is unclear whether old PNe like Sh 2-216 could be useful probes of the magnetized 
ISM. Targeted polarimetric observations at high resolution (<1') ; and possibly at multiple 
frequencies in the range 1-3 GHz, may help separate the signatures of more PNe from the 
turbulent ISM. 
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Table 1. Polarization Angles in the "Knots" of the North-East Arc 



Knot 


P.A. 


crp.A. 


Extrema 


n 


P.I. 


crp.i. 




(°) 


(°) 


(°) 




(K) 


(K) 




[1] 


[2] 


[3] 




[5] 


[6] 


1 


+33.6 


13.4 


+59.7, +21.0 


10 


0.139 


0.017 


2 


+83.7 


13.1 


-74.4, +66.0 


7 


0.125 


0.025 


3 


-88.1 


4.2 


-83.3, +86.8 


10 


0.122 


0.013 


4 


+52.7 


2.5 


+55.6, +48.1 


8 


0.229 


0.032 


5 


-79.9 


2.4 


-76.8, -84.3 


10 


0.175 


0.015 


6 


-76.6 


4.7 


-70.2, -82.5 


10 


0.217 


0.009 


7 


+89.9 


3.9 


-83.3, +84.9 


10 


0.187 


0.016 


8 


+65.7 


2.2 


+69.7, +62.0 


10 


0.220 


0.013 



Note. - - [1] Mean polarization angle over n pixels; [2] Stan- 
dard deviation in polarization angle over n pixels; [3] Low, high 
value of polarization angle (running clockwise from P.A. = 0); 
[4] Number of pixels used to estimate the mean polarization 
angle; [5] Polarized intensity at center of knot; [6] Standard de- 
viation in polarized intensity over n pixels. Polarization angles 
are modulo 180°; i.e., angles of —90° and +90° are equivalent. 
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159.5° 159.0° 158.5° 15B.0° 157.5° 159.5° 159.0° 158.5° 158.0° 157.5° 



GALACTIC LONGITUDE GALACTIC LONGITUDE 

Fig. I. — Images of the 2.5° x 2.5° region centered approximately (see text) on the position of 
PN Sh 2-216 in (a) optical intensity at R-band (from the DSS) and (ft) total radio intensity at 
1420 MHz. Here and hereafter, images are presented in Galactic coordinates, with Galactic 
north up and Galactic east to the left. The gray scale is in photon counts in (a) and brightness 
temperature in (ft), with lighter shades indicating higher counts/temperatures. The range 
of intensities in (a) has been adjusted to highlight extended emission. Point sources are 
removed in (ft) leaving extended emission with brightness temperatures in the range 4.57- 
4.92 K. The contours drawn in (b) accentuate higher brightness temperatures, and run from 
4.86 K to 4.92 K in steps of 0.01 K. The angular resolutions are (a) ~1" and (ft) 5' (smoothed 
from The dotted circle drawn on each image, and on each subsequent image, shows 

the approximate extent of the visible disk of Sh 2-216. 



-23- 





Fig. 2. — Images in (a) polarized intensity, P = >JQ 2 + U 2 — (1.2<r) 2 , and (b) polarization 
angle, 6p = \ arctan U/Q. The gray scale is in brightness temperature in (a) and runs from 
to 0.67 K, with lighter shades indicating higher temperatures. The gray scale in (b) extends 
from —90° (black) to +90° (white). Note that abrupt black-to- white transitions in (b) do not 
represent large changes in angle, since polarization angles of —90° and +90° are equivalent. 
The resolving beam in each image is 1.31' x 0.97' (full-width at half-maximum; FWHM) 
oriented at a position angle (east of north) of —40°. 
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Fig. 3. — Polarization angle image zoomed to a 0.4° x 0.4° region around the north-east 
arc. The gray scale is as described for Fig.©. The eight "knots" discussed in the text are 
labeled. The ellipses (thick dashed lines) drawn around each knot represent the resolving 
beam (FWHM), and define the knot perimeters for which the polarization angles (see Tabled]) 
are determined. 
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Fig. 4. — Simple model showing the interaction between Sh 2-216 and the ISM magnetic 
field. The perspective is that of an observer sitting above the Galactic plane and looking 
down on the center of the PN, with I = 158.5° directed upward. The PN is expanding in all 
directions. The motion of the PN system projected onto the Galactic plane is indicated. The 
host white dwarf ( "central" star) is offset from the center of the PN toward the interaction 
region, where the ISM magnetic field (solid lines) is compressed and deflected around the 
shell of the PN. Magnetic field lines which appear to stop at the edge of the PN actually 
slide on the surface either above or below the slice shown. The intrinsic ISM field is inclined 
15.5° to the plane of the sky. The lines-of-sight for an Earthbound observer to the north-east 
and north-west arcs are indicated. 



